function FreeBound(x,y)
if nargin==0
    global Node
    x=Node(:,2);
    y=Node(:,3);
end
dt=DelaunayTri(x,y);
ff=freeBoundary(dt);
hold on
plot(x(ff)',y(ff)','r')
